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Abstract 

We calculate the s— and d— wave superconductor order parameter in the 

vicinity of a tilt grain boundary. We do this self-consistently within the Bo- 

goliubov de Gennes equations, using a realistic microscopic model of the grain 

boundary. We present the first self-consistent calculations of supercurrent 

flows in such boundaries, obtaining the current-phase characteristics of grain 

boundaries in both s— wave and d— wave superconductors. 
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The debate over the superconducting order-parameter in the high T c superconductors 
(HTSC) has been strongly contested, but has now been settled in favour of a d— wave 
pairing state The experiments by Tsuei et al. and Wollman et al. [|J have been 

most conclusive, especially since they only depend on the phase of the order-parameter 
and not on the microscopic physics of the energy gap. Photoemission experiments []5| and 
the temperature depedence of the penetration depth || also strongly support the d— wave 
picture. However, some controversy still exists concerning possible s — d mixing in the 
cuprates. In particular an s-wave component has been demonstrated to be induced at 
interfaces 

In interpreting the Tsuei et al. experiments it is essential to understand the supercon- 
ducting characteristics of grain-boundary (GB) weak links in the cuprates. Understanding 
the effects of grain boundaries is also of importance for developing possible devices and other 
applications of high T c superconductors. The Tsuei et al. experiments, especially the ob- 
servations of 7r-j unction behaviour, were consistent with the predictions of d— wave pairing 
interpreted using the Sigrist-Rice []nj model for the dependence of the critical current, J c , 



on the grain boundary angles. On the other hand the values of I c measured as a function of 
grain boundary angle show an almost exponential dependence on angle [ [TT]JT^] , unlike the 
cosine predicted by the Sigrist-Rice formula for d— wave pairing. In trying to explain these 
findings a number of different models of the interfaces have been studied. Tanaka and col- 



legues []13 have looked at the (100) and (110) interfaces and also derived a Josephson current 
formula for s/I/d and d/I/d grain boundary structures. Barash et al. jl4| have considered 
the temperature dependence of the critical current in <i-wave junctions. Zhang considered 



0°, 45°, and 90° junctions [15]. Zhitomirsky and Walker [16] have also looked at the (110) 



interface to study the quasiparticle spectra and zero energy states (ZES). Beltzig et al. JT 
showed that an induced s— wave component existed on the (orthorhombic) (110) boundary 
giving rise to a splitting in the ZES at a low enough temperature: the latter point they 
attribute to Time Reversal Symmetry Breaking (TRSB). The review on GBs by Prester (TE| 
also highlights the possibility of them behaving as though each were an individual Joseph- 



son Junction. Gurevich and Pashitskii argued that the near exponential depencence of I c 
on angle was due to the formation of an insulating layer at the grain boundary, assosciated 



with the dislocation cores [19]. 

In this letter we address the effects of grain boundaries in both s— and d— wave supercon- 
ductors using a fully self-consistent solution of the Bogoliubov de Gennes (BdG) equations. 
We adopt a geometrically realistic model of the tilt grain boundary (GB), as shown in Fig.|l|. 
By solving the BdG equations in real space using the Recursion Method we are able to study 
such complex geometries, unlike earlier calculations which were limited to either simpler in- 



terfaces or planar junction models |PQ , |15| , pl]| ■ 111 consequence we can determine how the 
superconducting order-parameter (A), the charge density (n), and the quasi-particle local 
density of states (LDOS, n(E)) are affected by our GB. Further, by solving the BdG equa- 
tions in a self-consistent manner we can apply phase- differences in A across the boundary 
and calculate the resulting supercurrent. By calculating the maximum current across the 
boundary we determine the critical current of the system. 

For the purposes of this letter, we concentrate on the large angle grain-boundary (53.1°), 
depicted in Fig.l. It consists of two square lattices butted together at some angle of misori- 
entation and linked via a percolation site. Being periodic in the y direction we only have to 
undertake calculations for sites on two disimilar lines of atoms as shown in Fig. 1. For each 
line of sites we have to go 10 sites deep into the bulk before the order parameter has recov- 
ered to its bulk value. For some bonds across the GB the interatomic spacing is less than 
the bulk spacing. For these bonds we assume that the value for the hopping between sites i 
and j (tij) can be calculated by assuming a simple linear form for the hopping-integral, i.e. 

t* = --^ZT' - r * J - ^ (1) 

and is zero otherwise. For the geometry of our GB (Fig.l) every site will have a connectivity 
of 4. 

We consider the following two attractive Hubbard models: 
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with U < and n« = + n^. Here [/ is the usual BCS pairing-potential, defined as being 
a negative constant within a cutoff energy range of ±E C either side of the Fermi energy, 
after which its value is zero. Eq.2 will be refered to as 'local', giving rise to s— wave pairing 
and Eq.3 will be termed 'non-local' giving d— wave pairing. By making the Bogoliubov 
canonical- transformation we diagonalise the Hamiltonian and arrive at the Bogoliubov-de 
Gennes equations 
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where Ha = (— \i + \Uriu) and H i3 - = t i3 - (local), or, Ha = —jjl and = + \Uriij (non- 
local). Here \i is the chemical potential, it" and vf are the particle and hole amplitudes on 
site % associated with an eigenenergy E n , and is the appropriate charge density (defined 
below). To solve these equations we employ the Recursion Method [fffi] , and together with 
the methods employed in Martin and Annett [^3 , we obtain a matrix continued-fraction for 
the Green functions. This continued-fraction is evaluated exactly to 50 levels after which its 
elements vary in a predictable manner and therefore can be extrapolated for a futher 1500 
levels say. 

We are interested in evaluating the local quasi-particle density of states, the local and 
non-local charge densities {riu = J2a( c la c icr) an d n-ij = Do-( c !o- c j°-) respectively), and the 
local and non-local order-parameters (A& = U(c^Cii) and A^- = U(c^Cji) respectively). 
These quantities may be found from the Green functions, expressions for which have already 
been given elsewhere |[23|| . In the calculations that follow we make a BCS cutoff of U = 
— 3.5i, have a temperature, T, of T = O.Olt, and E c = 3.01 Iterating the equations for 
charge-densities and order-parameters, with the BdG equations, we generate self-consistent 
solutions. So as not to direct the final SC solution into a local energy minimum, we set the 



order-parameter to zero at the beginning of the calculation on those sites closest in proximity 
to the boundary. We say self-consitency has been achieved when the Hartree-Fock term and 
the order-parameter change by less than a predefined margin between iterations. For the 
s— wave case, we can reach 0.5% s.c. typically between 10 or 20 iterations, whereas for the 
d— wave it usually takes over 200 iterations. 

By imposing a phase- difference p in the order-parameter between the two bulk regions 
we may now generate current-flow across the GB. We initially make the Peierls substitution 
for tij, (Uj — > tije~ teAii *l r \ Aij being the integral of the vector-potential between sites i and 
j), in Eq.2 and 3, and use the definition 1^ = {■§¥-:)■, to obtain [^4] 
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f{E) being the Fermi-Dirac function. By a suitable choice of gauge we can immediately set 
Aij = everywhere. 

When applying phase-differences in A across a GB it is instructive to note how A changes 
due to the self-consitency, and also observe how the local densities of states and the Hartree- 
Fock term alter. Fig.^| shows the results for the s— wave case. Figs. 2(a) and 2(b) show the 
evolution of \Au\ across the GB on both the lower (a) and upper (b) lines of SC sites at 
<p = 0° (full line) and also (p = 180° (dotted line). The GB obviously has a perturbing effect 
on the system but note that |Ajj| is not depressed on the percolation site (x = 0) for p = 0°. 
Infact |Ajj| is almost constant as (p is varied; except for tp = 180° where |Ajj| = on the 
percolation site and is strongly reduced nearby. Note also the small Friedel-like oscillations 
in |Ajj| near the GB. Fig. 2(c) shows how the order parameter phase, arg(Au), varies as a 
function of x co-ordinate through the GB for the <p = 30° case. For the 180° case we just 
have a step function with arbitrary phase associated with the percolation site since |Ajj| = 0. 
For all other phase differences, the phase as a function of x goes as (y?/2) tanh(x/d) where 
<p is the bulk phase difference across the GB and d is a characteristic length scale which we 
find to be d = 3.5 for our parameters. Also, this equation holds for both SC lines. Finally, 



in Fig. 2(d) we show the local density of states on the percolation site at ip = 0° (full line) 
and at <p = 180° (dotted line). At 0° we have the usual BCS-like gap around the Fermi 
energy. This persists through all our calculated phase-differences except at 180° where the 
gap is suddenly filled in with various resonant states. Although not presented here we note 
that the Hartree-Fock term ^Unu shows the same (enhanced) values on both lines of SC 
sites at different phase-differences. 

Now consider the d— wave case. Because each site has four bonds with order parameters 
Ay we can calculate the 'net' d— wave contribution at a particular site by considering Af = 
^2j = i^(— l) J Ajj Similary, the extended— s component is given by A^ = Z)j=i,4 Ay and will be 
finite near the boundary because of the broken crystallographic symmetry. A minor difficulty 
now arises when discussing the d— wave scenario. In the definition of Af we must choose 
a convention for the direction of the positive and negative lobes of the d— wave function. 
For GB angles close to 45°, such as in Fig.l, there is an ambiguity in defining the relative 
orientations of the order parameter lobes either side of the GB. The convention we have 
chosen is indicated in Fig.l. In this case the current-phase relationship, 1(f), is similar for 
both s— wave and d— wave cases. For this choice the d— wave order parameter is zero on 
the percolation site, x = 0, for zero phase difference, Fig|5](a), unlike the s— wave case of 
Fig. 2 (a). Using our definition of phase-differences we find a maximum d— wave contribution 
at 180° on the lower line of SC sites (dotted line), which decreases with phase-difference, 
down to 0° (full line) where | A^ =0 | = 0. The upper line of SC sites (Fig.3(b)) shows the same 
qualitative form for A^ as for the local s— wave case. At 0° we find maximum extended— s 
component on the precolation site (full line in Fig.3(c)), decreasing with increasing phase- 
difference until at 180° we have minimum extended— s contribution (dotted line). Fig. 3(d) 
illustrates how the extended— s component evolves on the upper line of SC sites at either 
phase difference. Thus we conclude that the extended— s and the d— wave components are 
in competition such that the extended— s component is maximised at the detriment of the 
d— wave and vice versa. 

Fig.f| shows the calculated currents in both the local s— wave and non-local d— wave 
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cases. To calculate the current (using Eq.5) we have to consider the flow across all possible 
routes in just one cell of our sample. In calculating currents it is essential to check current 
conservation: this is only guaranteed from a self-consistent solution. Here, we find conser- 
vation obeyed to within 0.01%. Consider the s— wave intially: the variation of current with 
phase-difference, tp, is plotted in Fig. 4 (solid line) where the values for phase difference range 
from 0° to +180°. It is immediately clear that the current is not sinusoidal in ip but instead 
shows a sharp step at 180°. The remainder is roughly that of a saw-tooth albeit with some 
saturation. The step at ±180° can be attributed to resonant states entering the gap p5,26|. 



Fig. 2(d) compares the local quasi-particle density of states at <p = 0°cm<il80 o confirming the 
presence of resonant midgap states at 180°. 

Our calculations for the rf-wave current are presented in Fig.4 (dashed line). Again 
it is approximately a sawtooth. The slope |^ is also positive, and consequently this GB 
cannot be classified as a 7r-junction: this is consistent with the Sigrist-Rice formula for this 
geometry. Fogelstrom and Yip note that in certain geometries it is also possible to have 
a vanishing current at phase-differences other than integer multiples of tt, and this they 
attribute to time reversal symmetry breaking (this has been reported in Il'ichev P7[]). Fig.4 
shows no such evidence and therefore we conclude that the symmetric grain-boundary does 
not have TRSB. 

In conclusion we have developed a real-space method for determining how the order 
parameter and supercurrents change with phase-difference across a realistic interface in 
a superconductor. In this letter we have considered a large-angle symmetric tilt grain- 
boundary and considered the local s— wave and non-local d— wave pairing symmetry in the 
order-parameter on an equal footing. We have calculated the LDOS, Hartree-Fock terms, 
order parameter and current all self-consistently. By imposing a phase differnce, (p, across 
the junction we calculated the supercurrent We found, for both s— wave and d— wave 

that I(<p) is non-sinusoidal but exhibits a saw-tooth like behaviour which can be attributed 
to a sudden filling-in of the energy gap at tp = 180°. Further, we note no time reversal 
symmetry breaking or 7r-j unction behaviour in the d— wave case. 
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FIGURES 

FIG. 1. The symmetric model tilt grain-boundary. By periodicity, we carry out self- consistent 
calculations on two lines of sites (highlighted) which are then mirrored onto similar sites in the rest 
of the sample. Also shown is our definition of a d-wave 0° phase difference across the boundary 
(see text). 

FIG. 2. a) order-parameter on lower line of SC sites at 0° (full line) and 180° (dotted line). 
The latter goes to zero on the percolation site (x = 0), b)same as a) but for upper line of SC sites, 
c) self consistently determined evolution of the order-parameter phase on going through the GB for 
a 30° phase difference, d) Quasiparticle local density of states at 0° phase difference (full line) and 
at 180° (dashed line). 

FIG. 3. The contributions to the superconducting order-parameter (d—wave) on going through 
the boundary. The x coordinate is plotted on the horizontal axis. Full line refers to 0° phase 
difference, dashed refers to 180°: a) d-wave on lower line b) d-wave on upper line c) extended-s on 
lower line d) extended-s on upper line 

FIG. 4. The supercurent versus phase- difference in the order-parameter between the two bulk 
regions for the s-wave and d-wave cases. 
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